Momentum distribution and correlation function of quasicondensates in elongated 

traps 
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We calculate the spatial correlation function and momentum distribution of a phase-fluctuating, 
elongated three-dimensional condensate, in a trap and in free expansion. We take the inhomoge- 
neous density profile into account via a local density approximation. We find an almost Lorentzian 
momentum distribution, in stark contrast with a Heisenberg-limited Thomas-Fermi condensate. 
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Low-dimensional, degenerate Bose gases are expected 
to have significantly different coherence properties than 
their three-dimensional (3D) counterparts. In one- 
dimensional (ID) uniform systems, no true condensate 
can exist at any temperature T because of a large pop- 
ulation of low-lying states that destroys phase coherence 
(see |l| and references therein). For a trapped gas, the sit- 
uation is different: the finite size of the sample naturally 
introduces a low-momentum cutoff, and at sufficiently 
low temperature T <C X 1 ,/,, a phase coherent sample can 
exist lj. Above X</>. the degenerate cloud is a so-called 
quasicondensate: the density has the same smooth profile 
as a true condensate, but the phase fluctuates in space 
and time. As shown in , this analysis holds also for 3D 
condensates in elongated traps even if, strictly speaking, 
radial motion is not frozen. Such 3D, phase-fluctuating 
condensates have been recently observed experimentally 
in equilibrium [3j and nonequilibrium samples. 

Phase fluctuations of the condensate are caused mainly 
by long-wavelength (or low-energy) collective excitations 
IHu<lj| I n elongated traps, the lowest energy modes are 
ID excitations along the long axis of the trap 0. Fur- 
thermore, in the long-wavelength limit, density fluctua- 
tions are small and can be neglected for the calculation of 
the correlation function 0, 0] . Then, the single-particle 
density matrix is, assuming cylindrical symmetry, 



{&{p, zMp, z')) » X (P, Z, s)e-^^ 



(1) 



We have introduced A<j> 2 (Z,s) = ([<j>{z) - 0(z')} 2 ), the 
variance of the phase difference between two points z,z' 
on the axis of the trap, with mean coordinate Z = 
(z + z')/2 and relative distance s = z — z' , and the 
overlap function \ = \/ n o(p, z)n<j(p, z'), where hq is the 
(quasi)condensate density. The variance A<j) 2 (Z, s), the 
key quantity to characterize the spatial fluctuations of 
the phase of the condensate, has been calculated in 0, 
and an analytical form has been given, which is valid 
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near the center of the trap (i.e. for Z, s <C L, with L 
the condensate half-length). The first goal of this paper 
is to find an analytical approximation for the variance 
A<f) 2 (Z,s) valid across the whole sample. This is mo- 
tivated by the fact that coherence measurements with 
quasi-condensates El are quite sensitive to the inho- 
mogeneity of the sample. In position space, interferome- 
t rv [l3 . ll3j gives access to the spatial correlation function 
CW(s) (see e.g. |3) 

C«( S ) = J d 3 R(&(p, Z + s/2)*(p, Z - a/2)). (2) 

Equivalently, one can measure the axial {i.e. integrated 
over transverse momenta) momentum distribution V (p z ) , 
which is the Fourier transform of C^(s) 0HE[!J : 



V{Pz) = 2^h J dsC(1) ( s > 



-ip z s/h 



(3) 



A powerful tool to measure V(p z ) is Bragg spectroscoj. 
with large momentum transfer, as demonstrated in 
for a 3D condensate, and recently applied in our group to 
per fom the momentum spectroscopy of a quasicondensate 
[l0|. It is clear that both and V are sensitive to 
the inhomogeneity of the system. Our second goal is 
to obtain explicit expressions for these two important 
quantities. 

This paper is organized as follows. First, we summa- 
rize the results of [2| , and give an energetic interpretation 
of T^. Next, we discuss in detail a local density approach 
(LDA) to compute the variance of the phase for any mean 
position in the trap. This approximation is found to be 
accurate for T > 8X0, when applied to a trapped con- 
densate. Using the LDA, we then address the problem 
of a phase-fluctuating condensate in free expansion. In 
particular, we point out that at higher temperature, the 
phase fluctuations dominate over the mean-field release 
velocity and govern the shape of the momentum distri- 
bution. 

We consider A^o condensed atoms, trapped in a cylin- 
drically symmetric, harmonic trap with an aspect ratio 
A = lu z /u>± <C 1. If p ^> {hu>±,fru> z }, the condensate 
is in the 3D Thomas- Fermi (TF) regime 0. The den- 
sity has the well-known inverted parabola form: no(r) = 



2 



nom(l — p 2 — z 2 ), with the peak density nom = p/g re- 
lated to the chemical potential p. From now on, we use 
the reduced coordinates p = p/R and z = z/L, with 
R 2 = 2p/Moj 2 ± and L 2 = 2p/Muj 2 respectively. 

As shown in 1], phase fluctuations in trapped gases 
are mostly associated with thermally excited, low-energy 
quasi-particles (the quantum fluctuations are negligible). 
Under these conditions, the variance Acj) (Z, s) is 



Acj) (Z,s) 



where the sum extends over the ID axial excitations, with 
energy TkOj and occupation number Nj « k^Tjhaij for 
Nj 1. For a 3D condensate in an elongated trap, 
the amplitude cf>j is proportional to a Jacobi polynomial 

Pj , and ujj — LUz^jij + 3)/2 [f| for integer j. The 
explicit result for the variance is then Q 



Acj)\z,s) = ^ r f{Z,s), 



(5) 



with f(Z, S) = Fj [PW (Z + 5/2) - P^ l) (Z S/2)} 2 , 
and the coefficients Fj = (j + 2)(j + 3/2)/4j(j + 
l)(j + 3). Below the characteristic temperature T$ = 
15No(hu z ) 2 /32/i/cB, the phase profile is almost flat, and 
the single-particle density matrix is limited by the 
overlap function %: therefore the characteristic width 
of (i.e. the coherence length) is of order L. On 
the other hand, if T 3> 7$ the variance Acj> dominates 
the behavior of and the coherence length is sub- 

stantially smaller than L. Near the center of the trap 
(Z, s«l), Petrov et al. Q have derived the simple law 
Acj) (Z,s) ~ (T/Tj>)\s\, and introduced the characteristic 
phase coherence length = LT^/T that depends im- 
plicitly on the temperature, on the number of condensed 
atoms and on the trapping geometry. 

We can understand this expression for from ener- 
getic considerations. A random phase gradient of the 
condensate wavefunction, on a length scale L^, requires 
an average kinetic energy E<j, ~ N h 2 /ML 2 ^. This ki- 
netic energy is supplied by the thermal excitations that 
drive the fluctuations of the phase Q . As these excita- 
tions are quasi-classical (Nk 3> 1), this energy is of order 
k^T times the number of relevant modes. In ID /c-space, 
the distribution of the relevant excitations extends over 
~ and the spacing between modes is ~ l/L be- 

cause of the finite size of the system: this gives L/L^ 
relevant modes. By equating the two expressions for E^, 
we recover finally ~ LN (hu) z ) 2 / pksT . 

As indicated earlier, it is important to take the full 
spatial dependence of Acj) (Z, s) into account for quanti- 
tative comparison with experiments. In any case, Eq. @ 
can be evaluated numerically. However, we gain physical 
insight with an analytical approach based on the local 
density approximation (LDA), also used in [3j to calcu- 
late the evolution of the density in time-of- flight. This 



approximation considers that the condensate is locally 
equivalent to a homogeneous medium, however with a 
slowly varying density that depends on the trapping po- 
tential. If T ^S> T$, the coherence length is sufficiently 
small compared to L, that the LDA is valid for the cal- 
culation of correlation properties. 

The first step is to consider a finite cylinder of length 
2L, with radial harmonic trapping and periodic bound- 
ary conditions along z (and therefore homogeneous axial 
density). For this geometry, we find in the TF regime: 
no(r) = nom(l — P 2 ) for the condensate wavefunction. 
Low-lying excitations are found using standard Bogoli- 
ubov theory |5j, after averaging over the transverse de- 
grees of freedom [(|. The Bogoliubov spectrum for the 
excitation frequencies is — (D^i^k + 2-Mc^/fi.)) 1 / 2 f» 
cn)k for small k, with the free particle energy hiok : 
h 2 k 2 /2M and the ID speed of sound c iD = \J pj2M 
|19| . The Fourier component for phase fluctuations with 
wavevector k is 



2u k VV 



Mc 



ID 



1 



hk VV' 



(0) 



where the final expression holds for low-lying phonon 
states (k — > 0), and V = 2Tmo m R~ L. In a second step, 
we take into account the trapping potential by the sub- 
stitution: 



^MloIz 2 . 



(7) 



This implies directly the replacements: 



density : 
speed of sound 
radius : 
half — length : 



n m(oc p) 
ci D (cx y/p) 
R(oc y/Ji) 
L 



n0m(l 
ClD 

rVT 

L 



~z 2 ) 



(8) 



With these substitutions, we recover the 3D TF density 
profile. We require that the excitation frequency c\^>k 
is not modified as well, which implies replacing k with 
k(l — z 2 ) -1 / 2 , and using a density of states M{k)dk = 
(L/tt)(1 — z 2 )~ 1 / 2 dk. For the position-dependent vari- 
ance of the phase, we find [2(j 



(Z, s) 



T 



Tcf, (1 - Z 2 ) 2 



(9) 



The Z-dependent phase coherence length L$(l — Z 2 ) 2 
appearing in Eq. (jSJ can be substantially smaller near 
the edges of the trap than in the center. We will see that 
this reduces the average coherence length below its value 
at the center of the trap. 

We deduce from Eas. (|ll2l9|) the correlation function: 



T 



exp ( 



2T (1 - z 2 )- 



:)■ (10) 
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In deriving Eq. I|10|l , we have used the approximation for 
the overlap function x(p,Z, s ) ~ (1 — p 2 — z 2 — s 2 /4), 
valid near the center of the trap. In Fig. ^ we com- 
pare the result (|1(J|> to the correlation function following 
from the numerical integration of (JSJ. In the T = 
limit, the correlation function is limited by the overlap 
X- Because of the approximate form of Xi our result 
N Q (1 - (S/2) 2 ) 5 / 2 is about 25% too broad, 



C 



(i) 

trap,T=0 



and one should rather use the gaussian approximation 
to C^ 1 ' derived in 0. As T increases, C (1 ) turns to 
an exponential-like function, and our approximation ap- 
proaches the numerical calculation. For T > 8T^, the 
LDA result is very close to the numerical one (maximum 
error w 3%). For T > T$, Eq. JTUI) can be further sim- 
plified by keeping the s-dependent term only in the expo- 
nential. The Fourier transform then gives the momentum 
distribution 



V tla p,T(Pz) 



32tt 



dz- 



(l-Z 2 )4p2 + ^/4' 



(11) 



where p^, — h/L^ is a typical momentum associated with 
the phase fluctuations. This function is self-similar in 
Pz/p<t>, and approximated to better than 4% by a nor- 
malized Lorentzian with a half-width at half-maximum 
(HWHM) of Ap = O.67p0. This Lorentzian shape of 
the momentum distribution differs qualitatively from the 
fully coherent case, where it is almost Gaussian and lim- 
ited by the Heisenberg principle [14j . The increase of the 
phase fluctuations with increasing T not only broadens 
the momentum distribution, but also induces the appear- 
ance of " wings" , that form the "high-energy tail" of the 
quasicondensate. To quantify the accuracy of our ap- 
proximation, we have calculated numerically the Fourier 
transform of the correlation function. We find empiri- 
cally that the HWHM is accounted for by the formula 
Ap 2 w (2.04^/i) 2 + (O.65^/L ) 2 . The first term cor- 
responds to the Heisenberg-limitcd momentum width, 
and the second to the phase fluctuations. For T > 8T^, 
the height and width agree to better than 4% with the 
Lorentzian approximation. For lower T, the overlap func- 
tion x still affects the momentum distribution. 

Note finally that the momentum distribution is 
Lorentzian only in the domain k <C R . Outside of 
this region, the 3D nature of the excitations should to be 
taken into account properly. However, this does not af- 
fect the quasicondensate peak we are investigating here, 
but only the much smoother thermal background 

The results of the above paragraphs are valid for an 
equilibrium situation. However, coherence measurements 
involving Bragg scattering |l2l Il4| suffer from two ma- 
jor difficulties in a very elon gate d trap mean- field 
broadening of the resonance [l4|. and elastic scattering 
from the recoiling atoms and the condensate towards ini- 
tially empty modes fOj . Both of these problems can be 
solved by opening the trap abruptly, and letting the BEC 
expand to decrease its density before measurement. In 
the remainder of this paper, we discuss how expansion 
modifies the momentum distribution and the correlation 



function, assuming that the expansion time is chosen to 
be long enough to neglect the collisions. 

For a pure, elongated condensate abruptly released 
from the trap at t = 0, the explicit solution was found 
in 22] . The condensate density keeps its initial Thomas- 
Fermi shape, with the coordinates re-scaled. The (small) 
axial momentum from the released mean-field energy is 
linear in position: p z w p cxv z, with p cyLp — (ir /\/2)\Mcs 
for t — oj±t 1, and eg — y/Ji/M is the 3D speed 
of sound. The axial momentum distribution mirrors the 
(integrated) density distribution: 



V, 



cxp.T: 



15 



16p, 



OXp 



, Pz - 
Pcxp 



(12) 



This expression holds for a pure condensate, at T = 0, 
as indicated. For a phase fluctuating condensate at finite 
T, it is necessary to consider the time evolution of the 
fluctuations as well. As shown in 3], the momentum 
distribution partially converts into density modulations 
after time of flight. An explicit solution was derived for 
the density fluctuations in the axially homogeneous case. 
Using the continuity equation (after radial averaging) , we 
find for t 3> 1 



<pk(z,r) « cf>k(z, 0)t 



-(a 



^) 2 cos (^ r) , 



(13) 



If u z t < (/i/7kj z )(T/T ) 2 , then for all k < L^ 1 , the phase 
distribution is essentially frozen: <j)y{z,T) w 0k(z,O). 
Physically, this condition states that for such a time of 
flight, the excitations that have significant contributions 
to the phase fluctuations have not yet been converted into 
density modulations. This condition is not at all restric- 
tive for typical experimental parameters 0, 0, IToL Hl| . 
and we suppose it is met in the remainder of the paper. 

Using the rescaled wavefunction from (23, together 
with ®, we find the correlation function for the expand- 
ing quasicondensate: 



7T/1 T I S I 

eXp( ^ S -27>(I^)- ( 14 ) 

The phase factor in (|14f> accounts for the local expansion 
momentum introduced above (recall nfi/huj± 3> 1). The 
Fourier transform gives the momentum distribution 

7>oxp,T(Pz) « ^-ff 7=p ^/ Poxp (^). (15) 

Pcxp Pcxp 

The function g 1 is given by 



32tt7„ 1 (1 _ f2)4( x _ ~ z y + YL 



(16) 



and the ratio 7 = p<p/p cxp controls which component 
of the momentum distribution dominates. In the limit 



4 



T — > 0, using 7/(x 2 + 7 2 ) — > ttS(x) as 7 — > 0, we re- 
cover the zero-temperature result Eq. (|12fl. On the other 
hand, if 7 3> 1, we expect the momentum distribution 
to be similar to the distribution in the trap iJUJ. Fig. [21 
shows a numerical calculation of <? 7 , for various values 
of 7. We find that already for p^ > 2p oxp , the momen- 
tum distribution is almost entirely dominated by phase 
fluctuations, and, as in the trapped case, is very well 
approximated by a normalized Lorentzian with HWHM 
= 0.67^. Here, we note two points: first, that the 
Heisenberg width ~ %/L is negligible at any tempera- 
ture, and second that, for large enough condensates, we 
can have p cxp 3> p^, even if the coherence length is smaller 
than L. 

In summary, we have analyzed the measurement 
of phase fluctuations in elongated Bose condensates. 
Within a local density approach, we have been able to 
take the density profile into account, and derived analyt- 



ical formulas for the correlation function and the momen- 
tum distribution of static and freely expanding quasicon- 
densates. In the regime of interest, the formula compare 
well to a numerical evaluation based on the results of , 
which are exact in the long- wavelength limit. In partic- 
ular, we show how the shape of the momentum distribu- 
tion tends to a Lorentzian with half-width ss 0.67h/L,p 
as one goes further in the phase-fluctuating regime. We 
believe that these results may be helpful to understand 
quantitatively experiments involving quasi-condensates 
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FIG. 1: Spatial correlation function C^'(S) of a trapped qua- 
sicondensate. The solid lines follow from the numerical eval- 
uation of the result © derived in @|, for T = 0, T/T^ = 4 
and T/T,), — 8, in order of decreasing width. The dashed lines 
follow from the local density approximation Ea. l|10^ . Note 
that for convenience, C^(s) as shown here is normalized to 1 
rather than to No, differently from what is done in the text. 
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FIG. 2: The momentum distribution in expansion g 7 (Eq. ll6H 
for 7 = p^/pcxp = 0,1,3. As 7 increases, (? 7 continuously 
transforms from a quartic profile to a lorentzian-like profile 
(see text). For 7 = 3, we find little change from the momen- 
tum distribution in the trap (dashed line, a Lorentzian with 
HWHM 0.677, see text). The functions have been rescaled 
by their maximum values to facilitate comparison. 



